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Abstract. Photomeson production is the main energy loss for relativistic nuclcons 
in dense radiation fields like the cosmic microwave background and the radiation 
fields in Gamma Ray Bursts (GRB) and jets of Active Galactic Nuclei (AGN). In 
this paper we study photomeson production in typical GRB and AGN jet radiation 
fields by using the recently developed Monte Carlo event generator SOPHIA (see 
these proceedings). We discuss processes that are relevant for the physics of cosmic 
ray acceleration and the production of neutrinos and gamma rays. We compare our 
results with widely used approximations, and find significant deviations, particularly 
for GRBs. The photoproduction of antibaryons as a so far not considered effect in 
astrophysics is briefly discussed. 



1. Introduction 

Ultrahigh-energy nucleons propagating through dense radiation fields lose their 
energy mainly through photomeson production. Radiation fields are universal 
(Cosmic Microwave Background Radiation (CMBR)) and are very intense at luminous 
astrophysical sources such as Active Galactic Nuclei (AGN) jets and Gamma Ray 
Bursts. The secondary products of photomeson interactions decay, and lead eventually 
to the emission of neutrino and gamma-ray fluxes from the source, which may be 
observable. This has triggered several authors (e.g. Stecker 1973 ||, Berezinsky & 
Gazizov 1993 [|J, Mannheim & Schlickeiser 1994 j|) to study photopion production 
in more detail. Collider measurements indicate a complicated structure of the cross 
section, especially in the astrophysically important lower energy range. In order to 
overcome this problem, either simplified approximations for the cross section, e.g. the 
so-called A-approximation (Stecker 1973 [||, Gaisser et al. 1995 Q), or Stecker's 
isobar model have been used. Berezinsky & Gazizov || used interpolated cross 
section measurements and implemented them into numerical codes to derive neutrino 
production spectra. Recently, it has been shown that in realistic photon fields photon- 
proton collisions may happen at low and high center-of-mass (CM) energies (Miickc 
et al 1999a ||), making a more realistic treatment of the pion production process 
necessary. For example, hadronic interactions of ultrarelativistic protons in flat 
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ambient photon spectra, e.g. in Gamma Ray Burst radiation fields and the synchrotron 
radiation field in 7-ray loud AGN, are roughly equally likely to occur in the resonance 
and the multiparticle production regions. 

In this paper we discuss photopion production in typical GRB and AGN jet 
radiation fields. We discuss microphysical quantities relevant for the photohadronic 
production of gamma-rays, neutrinos, and neutrons (Section 2), by using our 
newly developed Monte Carlo event generator SOPHIA (presented also in these 
proceedings (Miicke et al 1999b 0)), and compare them with results from a widely 
used approximation (Section 3). As a qualitatively new effect we also discuss the 
photoproduction of anti-baryons (Section 4). The relevance of our results is briefly 
discussed in the context of neutrino and cosmic ray production (Section 5). 



2. Relevant microphysical quantities 

One of the major motivations for discussing photohadronic interactions in 
astrophysical sources like AGN jets and GRB is the prediction of observable neutrino 
fluxes from these systems. The total flux of neutrinos for a specific source model 
can be related to the predicted fluxes of gamma rays and cosmic rays, which in turn 
can be compared with current observations. Neutrino fluxes from hadronic AGN 
jet models are not expected to be observable as point sources by current neutrino 
telescopes, and therefore contribute to the extragalactic neutrino background. Hence, 
estimates for diffuse neutrino fluxes originating from blazar jets can be derived from 
measurements of the diffuse extragalactic gamma ray background (EGRB) (see e.g. 
Mannheim 1995 [||) if the entire EGRB is produced photohadronically. Upper limits 
can be set if only part of the EGRB is due to unresolved blazars (see Chiang & 
Mukherjee 1997 §, Miicke & Pohl 1998 0). 

Of particular relevance is here the production of gamma rays. In photohadronic 
sources, the opacity for the primary gamma rays emerging from the decay of neutral 
pions must be large, since otherwise the efficiency for the photohadronic interactions 
themselves would be too low to produce observable fluxes. The interactions of primary 
photons (and electrons) cause electromagnetic cascades which reprocess the leptonic 
power to lower photon energies, and photons are eventually emitted in the range < 100 
GeV (Mannheim 1993 (n)). 

The exact relation between the gamma ray and neutrino fluxes is determined by the 
microphysics of photomeson production, i.e., by the photon-to-neutrino total energy 
ratio per interaction, S 1 /E v — Y^E^/Y^Ey. The A-approximation, which has been 
often used, gives £ 7 /£„ = 3. This ratio has been slightly modified by considering 
other photon producing processes, like Bcthc- Heftier pair production off the energetic 



protons (Mannheim 1993 |11 ) 



Another way to constrain neutrino fluxes is a comparison with the cosmic ray 
spectrum. This is motivated by the fact that neutrino producing interactions also 
cause an isospin-flip of the proton into a neutron, which is not magnetically confined 
and can be ejected from the accelerator. In contrast, protons need to be confined in 
order to be accelerated, and their escape from sources involving relativistic flows, like 
AGN jets or GRB, is strongly affected by adiabatic losses. Neutron ejection can also 
be suppressed if the opacity for interactions exceeds unity. It is possible however to 
relate the 717 opacity to the 77 opacity of very high energy photons, and one can show 
that sources transparent to very high energy (VHE) gamma rays are also transparent 
to neutrons at most energies [O. Then, regardless of the proton confinement, the 
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minimum cosmic ray ejection is given by the produced neutron flux, which allows an 
upper limit on the possible neutrino flux (Waxman & Bahcall f999 It has been 

pointed out by Mannheim et al. ]l3fl , that the complicated propagation properties of 
cosmic rays make it difficult to apply a model independent cosmic ray bound on the 
neutrino flux. The relation of the cosmic ray flux to neutrino flux is again determined 
by a microphysical parameter, namely the neutrino-to-neutron total energy relation 
£v/£n = E n , which is predicted ~ 0.19 in the A-approximation. 

Other kinematical quantities of interest are the inelasticity of the proton, AE p / E p , 
which determines the proton energy loss time scale due to photohadronic interactions, 
and the average fractional neutrino energy per interaction, (E v ) /E p , which determines 
the maximum neutrino energy if the maximum proton energy is determined by the 
acceleration process. In the same way, we may also consider the fractional neutron 
energy per interaction. 

3. Cosmic ray, gamma ray and neutrino production 

In the typical scenario of hadronic interactions in AGN jet and GRB radiation fields, 
relativistic nucleons with energy E p 3> e are assumed to interact in an isotropic 
radiation field in the comoving frame of the relativistic plasma flow. Because of 
the strong magnetic fields required to accelerate protons, photoproduced electrons 
and positrons can be assumed to have 100% radiation efficiency, and transform their 
energy completely into photons through synchrotron emission. The total and average 
energy of the produced neutrinos includes v e and and their antiparticles. The 
v^jv^ contribution, possibly detectable in current or planned underwater/-ice neutrino 
experiments, is roughly 2/3 of the total neutrino flux. We assume that protons are 
magnetically confined, and count only photohadronically produced neutrons as ejected 
cosmic rays. 

3.1. Blazar/AGN jets 

There exists a large variety of models to explain the observed 7-ray emission from 
blazars (= flat spectrum radio quasars and BL Lac objects). The leptonic models, 
which assume inverse Compton scattering of low energy photons up to gamma ray 
energies, currently dominate the thinking of the scientific community. Alternatively, 
photopion production has been proposed to be the origin for the observed 7-ray flux 
(Mannheim 1993 jll]). A clear distinction between the two models is the production 
of neutrinos, which is negligible in the leptonic models, but may be equally important 
as the photon production in the hadronic models, and has been predicted to cause 
observable fluxes in large underwater/-ice neutrino detectors. 

Gamma ray and neutrino production in hadronic blazar models occurs through 
photopion production (and subsequent cascading) of relativistic protons in cither 
external (accretion disk or the IR radiation field from a molecular torus) radiation 
fields or the synchrotron radiation field produced by electrons which arc shock 
accelerated together with the protons. Here, we confine our discussion to the latter 
model. For simplicity, we consider TeV-blazars only, which are characterized by a low 
energy spectrum extending to X-ray energies. The synchrotron emission from blazar 
jets at low energies can be well explained by a superposition of several self-absorbed 
synchrotron components (e.g. Cotton et al 1980 Q), Shaffer & Marscher 1979 
leading to a flat target spectrum. The local synchrotron radiation spectrum follows 
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a power law with photon index a = 1.5 (n(e) cx e~~ a ) up to a break energy (see e.g. 
Mannheim 1993 pi). Above the break energy ef, « 10~ 4 eV (see e.g. Rachen, & 
Meszaros pi) it is loss dominated to become n(e) oc e~ 2 . The synchrotron emission 
in TeV-blazars like Mkn 421 and Mkn 501 is observed to continue up to 1-10 keV. 
Thus, for our application we approximate the typical target photon spectrum in the 
jet frame (assuming a Doppler factor of D=10) by 

n(e) cx e~ 3/2 for 10~ 5 eV < e < 10~ 4 eV 

n{e) cx e~ 2 for 10~ 4 eV < e < 10 3 eV (1) 

For the maximum proton energy we use -E^max = 7p«ipC 2 = 10 10 GeV in the rest 
frame of the jet, corresponding roughly to the highest observed cosmic ray energies 
after Doppler boosting. This is also consistent in order of magnitude with the values 
derived by equating the acceleration time of the proton with its total loss time due to 
adiabatic, synchrotron and photohadronic cooling (see Rachen & Meszaros |l6j for a 
more detailed discussion). Due to the threshold condition for photopion production 
of protons with Lorentz factor j p interacting with photons of energy e 

only photons from the steep part of the target radiation held interact photohadroni- 
cally with protons of 7 P < 10 10 . 

Fig. |j] shows £ 1 j£ v from SOPHIA simulations for protons of energy E p interacting 
in the synchrotron field of blazar jets. The prominent resonance/threshold region 
dominates the interaction, independently of input proton energy. Consequently, the 
total fractional photon energy £j/E p sa 0.1, independent of E p , and an overall photon- 
to-neutrino ratio of £~/£ v ~ 1.2. We note that this result differs by a factor ss 3 from 
the prediction of the A-approximation, although the steep spectrum emphasizes the 
threshold region with the dominant A(1232) resonance. This can be understood in 
view of the strong contribution of non-resonant direct ir + production at threshold, 
and the contribution of other resonances with different isospin in the vicinity of the 
A(1232) resonance (see Miicke et al., these proceedings [Q). 

In contrast, kincmatical quantities like inelasticity and fractional secondary 
particle energy are well reproduced by the A-approximation (see Fig. ||), since these 
quantities vary only slowly with energy. For example, the proton inelasticity is in 
general logarithmically rising with increasing nucleon input energy in the resonance 
region (Miicke et al 1999b Q). Thus, when interactions occur near the A(1232)- 
resonance like in this case, the A-approximation gives a good description. The same 
applies for the total neutrino-to-neutron energy ratio. Here, we find £ v j£ n ~ 0.18 (see 
Fig- 1). 

3.2. Gamma Ray Bursts 

The cosmological Gamma Ray Burst fireball model has been very successful in 
explaining the observed temporal evolution of the afterglow photon spectra (see 
Piran 1998 jl7| and references therein). In this model a relativistically expanding 
fireball transforms most of the explosion energy into the kinetic energy of baryons 
in a relativistic blast wave with bulk Lorentz factor of 100 - 300. This energy is 
reconverted into radiation at shocks, produced either in collisions between different 
shells ejected from the central source (internal shock scenario), or through deceleration 



40 



■5 

4 




Figure 1. The total photon and f-power (normalized to the proton input energy 
E p ), and their ratio emerging from a source due to photopion production of protons 
of energy E p in a typical TeV-blazar ambient radiation field The E 1 IS V — ratio is 
approximately unity, independent of proton input energy, which can be contrasted to 
the A-approximation (dotted lines). Solid lines denote the results from the SOPHIA 
simulations. 



of the blast wave when it sweeps up the external medium (external shock scenario). 
The former process is thought to be responsible for the GRB itself, while the latter 
produces the afterglow. In this scenario, the observed photon radiation is explained 
as synchrotron radiation from the accelerated electrons. In the comoving shock frame 
this target radiation field for photopion production of the relativistic protons may be 
approximated by a broken power law: 

n(e) cx e~ 2/3 for l(T 3 eV < e < 10 3 eV 

n(e) cx e- 2 for 10 3 eV < e < 10 5 eV (3) 

Waxman |ll| and Vietri |l9| suggested that if a significant fraction of the observed 
GRB power is transformed into ultrahigh energy cosmic rays (UHECRs), GRB may 
well be the source for all the observed UHECRs. In fact, it has been shown that 
comoving proton Lorentz factors of up to ~ 10 9 can be reached in GRB shells, which 
are boosted in the observers frame to 3 • 10 20 eV (||T^],p6|). This scenario gives 
also rise to fluxes of very high energy neutrinos (E v > 100 TeV) correlated with 
gamma ray bursts, which could be detected in a km 3 neutrino observatory because 
the exact temporal and directional information reduces the background to virtually 
zero (Waxman & Bahcall 1997 pTj[|). The neutrino emission, however, is suppressed 
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Figure 2. Average neutron (upper curves) and neutrino energy (lower curves) 
with respect to the proton input energy E p , and the mean proton inelasticity due to 
photon-proton pion production of protons with energy E p in a typical TeV-blazar 
ambient radiation field. The dotted lines represent the respective values from the 
A-approximation while solid lines denote the results from the SOPHIA simulations. 



at the highest energies because of the adiabatic and synchrotron losses of pions and 
muons prior to their decay (Waxman & Bahcall 1997 pOfl , Rachen & Meszaros 1998 
@, Waxman & Bahcall 1999 @). 

In the highly energetic GRB photon field mostly photons from the fiat part of 
the ambient radiation field interact photohadronically. For incident protons with 
7 P > 10 7 interactions predominantly occur in the multiparticle region of the cross 
section, while for 7 P < 10 7 the resonance/threshold region determines the secondary 
particle production. This leads to an increase of photon and neutrino production by 
roughly a factor of 2 with £ 7 w for protons with -f p > 10 7 in comparison with 
protons with j p w 10 5 (see Fig. |J). The photon-to- neutrino total energy ratio seems 
to be fairly robust £ 7 w £ v k, 1 at all proton energies, except for a slight deviation at 
7 P « 3 • 10 5 where photons from the flat part of the target photon spectrum interact 
mainly via the A(1232)-resonance. Also here, this result is about a factor of 3 different 
from the A-approximation. 

The average number of neutrinos produced per interaction can increase with CM 
energy by up to an order of magnitude for the proton energy range relevant in GRB (see 
Miicke et al 1999b 0). Because of the high multiplicity of the secondaries, the mean 
energy of the neutrinos produced in the multipion region (i.e. for j p > 10 7 ) is « 1% 
of the proton input energy (see Fig. |J). At lower energies it reaches (E„)/E p w 0.04, 
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Figure 3. The ratio of total ^-energy to neutron energy emerging from a TeV-blazar 
jet due to photopion production of a protons with energy E p . The A-approximation 
(dotted line) gives an adequate description of the v-to-neutron ratio. 



approximately in agreement with the A-approximation. Therefore, the ejection of 
ultra-high energy neutrinos (E v > 10 19 eV), as recently proposed by Vietri pl[ 
(but note also |l2| ) requires observer frame proton energies of 10 21 eV. Such high 
proton energies are not expected from shock acceleration scenarios in GRB shells. In 
the general case, however, the synchrotron losses of charged pions and muons in the 
highly magnetic (B ~ 10 3 — 10 6 G) GRB environment and their adiabatic deceleration 
determine the high energy end of the observable neutrino flux from GRB (see Rachen 
& Meszaros 1998 

If protons are magnetically confined at the source, the proton-neutron conversion 
probability determines the production of cosmic rays in photohadronic sources. At 
high energies it decreases from w 0.5 to 0.3 (see Miicke et al, these proceedings 0), 
and the average fractional energy of neutrons decreases to < 0.5 of the proton energy. 
Together with the increasing power going into secondary particles, this leads to a 
neutrino-to-neutron ratio £„/£ n ~ 1, about a factor of 5 above the A-approximation. 
This increase of the neutrino power relative to the cosmic ray power of GRB (if protons 
are confined) is, however, masked by the secondary particle losses (see Rachen & 
Meszaros 1998 [jf6| , Waxman & Bahcall 1999 [fl2|| ). At lower energies the neutrino-to- 
neutron ratio approaches the value expected in the A-approximation (see Fig. ||). 
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Figure 4. The total power in photons and neutrinos (normalized to the proton 
input energy E p ), and their ratio due to photomcson production in a typical GRB 
ambient radiation field. The £ 7 /£„— ratio is nearly always unity, independent of 
proton input energy, which can be contrasted to the A-approximation (dotted lines). 
Solid lines denote the results from the SOPHIA simulations. 



4. High energy anti-baryon production 

A qualitatively new effect which can be investigated with SOPHIA is the production of 
secondary baryon/anti-baryon pairs in high energy photohadronic interactions. This 
requires a detailed simulation of QCD string hadronization, which is implemented 
in SOPHIA. The theoretical threshold of this process is y/s — 3m p c 2 , corresponding 
to a photon energy in the proton rest frame e' = ArripC 2 . For e' > 2 TeV, about 
40% of all events contain antinuclcons. This raises the interesting question about 
the contribution of antiprotons to the cosmic ray flux from extragalactic photopion 
production sources. 

Due to the threshold condition, this process can only affect the antimatter to 
matter ratio at extremely high energy. Assuming the minimal cosmic ray ejection 
hypothesis, (only neutral baryons can leave the source and become cosmic rays) the 
relevant quantities are the anti- neutron to neutron ratio £n/£ n and the corresponding 
multiplicity ratio Nn/N n . 

As noted above, cosmic rays with energy 7 P > 10 7 accelerated in GRBs tend 
to interact with photons from the flat part of the ambient radiation fields, in the 
multiparticle production regime. While the total energy dissipated into photons, 
neutrinos and neutrons is roughly the same (£ 7 m £ v w £ n w Q.2E P ) in this high 
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Figure 5. Average neutron and neutrino energy with respect to the proton input 
energy E p , and the mean proton inelasticity due to photon-proton pion production 
in a typical GRB radiation field. The dotted lines represent the respective values 
from the A-approximation which are only met for interactions near threshold, e.g. 
at low input energies 7 P < 10 6 . Solid lines denote the results from the SOPHIA 
simulations. 



energy part of the cross section, the antineutron production is only 1/40 of the neutron 
production. The neutrons carry on average 50% of the input energy (see Fig. ||) with 
a neutron multiplicity of roughly 0.4. Antineutrons have mean energies of ~ 0.1E p 
with a multiplicity of approximately 0.05. Antineutrons and neutrons of the same 
energy are produced by protons of different energy. For example, we need an input 
proton energy of 10 10 GeV for the production of a antineutron of energy 10 9 GeV while 
for a neutron of the same energy the proton input energy must be ~ 2 • 10 9 GeV. One 
has to weight the ratio of the resulting neutron and antineutron multiplicities, and 
this is determined by the input proton spectrum. For this purpose we use a E~' 2 
differential spectrum as a typical equilibrium particle spectrum for GRB. For photon- 
proton interactions at high CM energies, as is the case for GRBs, we therefore expect 
antineutron-to-neutron multiplicities of the order 0.01. 

Since photon-proton interactions in TeV-blazar ambient photon fields take part 
predominantly at threshold, these objects are not expected to be strong sources 
of antinucleons. SOPHIA simulations show that the total energy for antineutron 
production does not exceed 10~ 4 of the proton input energy with the average 
antineutron carrying roughly 1/10 of E p . Neutrons, on the other hand, take up about 
20% of the input energy with a mean particle energy of 0.8E p , as expected from the 



10 

1 0.0 



LU 




0.1 



1 



5 



1 6 10 
Log(7 



7 R 9 10 

/ 10 ° 10 



e /D ) 

P 300' 300 y 



Figure 6. The total v-energy to neutron energy emerging from a source due to 
photopion production in a typical GRB photon field. With increasing input proton 
energy the total neutrino power per produced cosmic ray is underestimated by a 
factor < 5 in comparison to the A— approximation (dotted line). 



A-approximation. Weighting again with the typical equilibrium differential particle 
spectrum expected in blazar jet environments, which follows a E^ 1 power law, we find 
that the fi/n-ratio does not exceed 10~ 3 . 

The n/n-ratios resulting from SOPHIA in GRB and TeV-blazar radiation fields 
as a function of proton energy are shown in Fig. 0. After /3-decay, the ejected n/n 
contribute to the overall p/p flux. The n/n ratio derived here has to be regarded as an 
upper limit to the observable p/p ratio, since direct proton injection may contribute 
to the cosmic ray flux without increasing the anti-proton population. 



5. Discussion and outlook 



Our simulation results obtained with the new photomcson production event generator 
SOPHIA jjj demonstrate that the widely used A-approximation for photohadronic 
interactions has only a limited applicability to many important astrophysical 
applications, namely the secondary gamma ray, neutrino and neutron production in 
radiation fields of AGN jets and GRB. For AGN jets, where interactions at low center- 
of-mass energy dominate, the major kinematical quantities, like proton inelasticity, 
average fractional neutrino energy per interaction, and the neutrino-to-neutron energy 
ratio are found to be well predicted by the A-approximation. For interactions of ultra- 



11 




Figure 7. Predicted antineutron-to-neutron ratio at production in the shock frame 
of Gamma Ray Burst (GRB) and in TeV-blazar jets. GRBs predict a antibaryon- 
to-baryon production ratio, which approaches 0.01 above 10 eV, while TeV-blazars 
may produce antibaryons on a level of about one order of magnitude lower. 



high energy cosmic rays in GRB shells, however, we find deviations up to an order of 
magnitude for some of these quantities. 

As a robust result for both kinds of sources we find that the average total energy 
channeled into neutrinos, £ Vl and gamma-rays, £ 7 are approximately equal. By 
comparison, the A-approximation predicts £ 7 /£„ ~ 3. This relation has been widely 
used to normalize expected neutrino fluxes from AGN jets to their photon flux. As a 
consequence of this deviation from the A-approximation found in our SOPHIA results, 
the expected neutrino fluxes from such AGN models would increase by a factor of 3 
assuming that the observed gamma ray emission is entirely of photohadronic origin. 
The lower gamma-ray-to-neutrino ratio seen in our SOPHIA results, together with the 
neutrino-to-neutron ratio in agreement with the A-approximation, implies also that 
AGN neutrino models, which were initially designed to comply with the limitation 
set by the cosmic ray data in a straight line propagation scenario (Mannheim, 1995 
||, Model A), can no longer produce a large fraction of the 100 MeV gamma ray 
background as initially assumed. This may comply with the recent finding of Chiang & 
Mukherjee (1997 ||) that only part of the measured EGRB may be due to unresolved 
blazars. However, it has been pointed out by Mannheim et al. |l3| that an increased 
and so far unexpected contribution from sources with maximum cosmic ray energies 
below 10 19 eV can evade the problem without being in conflict with the cosmic ray 
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data. Moreover, it is possible that the optical depth of blazars for ultra-high energy 
neutrons is larger than unity, which allows a higher diffuse neutrino and gamma ray 
flux for a given cosmic ray flux. 

For the derivation of the upper bound of the neutrino emission, the neutrino-to- 
neutron energy ratio reflects the microphysics of the photohadronic interactions. The 
value EujEn ~ 0.19, expected from the A-approximation and used in Waxman and 
Bahcall [Q, was confirmed for AGN jets in our simulations. For GRB, however, 
we find a ratio of about 1. This implies that the cosmic ray upper bound for GRB 
neutrinos has to be increased by a factor of 5 at high energies. In practice this 
theoretical bound is unlikely to be reached by GRB neutrinos, since their emission 
is suppressed at high energies by secondary particle (pion and muon) cooling in the 
hadronic cascade [l6| . Models which try to evade such losses by shifting the 
acceleration region into the outer shock and the afterglow of the GRB, as suggested 
by Vietri ( PH ; note also Jl^|), can be shown not to be efficient enough to reach such 
fluxes for reasonable energetics |fT6f . 

The average neutrino energy due to photopion production in GRB spectra at very 
high proton energies is up to an order of magnitude below the value expected from the 
A-approximation. This has severe consequences for the prospects of GRB-correlated 
neutrino events above 10 19 eV. The fact that the mean proton energy loss rises with 
increasing proton energy up to AE/E ~ 0.5 (compared to the A-approximation 
value AE/E = 0.2) may add to the problem, since it leads to lower proton energies 
in photohadronically limited acceleration scenarios. 

A qualitatively new feature which can be explored with SOPHIA is the 
photoproduction of antibaryons |Q. Our prediction of the maximal antiproton 
contribution from GRB of about 2%, reached for energies above 5 • 10 17 eV, is several 
orders of magnitude above the background expected at this energy from cosmic ray- 
nucleon collisions in the galactic disk. The latter is measured at energies ~ 10 GeV 
as approximately 10~ 3 , and expected to decrease with E p '°' 3 °' 6 ' following the leaky 
box model. Since the expectation of a significant increase of the p contribution at 
high energies is unique for GRB, this could provide an independent test whether GRB 
are indeed the dominant sources of cosmic rays. Unfortunately, there is presently 
no imaginable way to distinguish between nucleons and anti-nucleons at such high 
energies, where cosmic rays can only be measured in air shower experiments. Our 
result may therefore be regarded as of rather academic interest, or may be kept in 
mind for presently not anticipated, future detection techniques. 

In conclusion, we have demonstrated that the application of SOPHIA to 
astrophysical problems involving the interaction of energetic cosmic rays in photon 
backgrounds can (a) improve the accuracy of the predictions from such models, and (b) 
open the possibility to explore particle physics effects so far neglected in astrophysics. 
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